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Abstract 

We compute the form factors of exponential operators e k9 ‘ p< ' x ^ in the two-dimensional integrable Bullough 

to\ 

Dodd model {a\ J Affine Toda Field Theory). These form factors are selected among the solutions of general 
nonderivative scalar operators by their asymptotic cluster property. Through analitical continuation to 
complex values of the coupling constant these solutions permit to compute the form factors of scaling 
relevant primary fields in the lightest-breather sector of integrable <j> 1,2 and cj) 1,5 deformations of conformal 
minimal models. We also obtain the exact wave-function renormalization constant Z(g) of the model and 
the properly normalized form factors of the operators tp{x) and : ip 2 ( x ): . 
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1 Introduction 


In recent years much important progress has been achieved in the study of two dimensional Quantum Field 
Theory and related Statistical Mechanical systems. The solution of Conformal Field Theories [[l| || has not only 
allowed the full characterizaton of fixed points in the renormalization group describing universality classes of 
critical models, but it has also provided the possibility of describing the renormalization group flow away from 
criticality by means of relevant deformations of Conformal Minimal Models |]]. In particular, Zamolodchikov 
showed that in some interesting cases Q falling in the classes of 0i,3 and 02,1 deformations — in¬ 
finitely many integrals of motion survive the deformation and the system is suitably described by an integrable 
relativistic scattering theory. Bootstrap techniques relying on the integrability of the model provide then a 
powerful tool for obtaining the exact 5-matrix of the system which turns out to be elastic and factorizable 
&§■ A systematic description of the 5-matrices for all the integrable deformed minimal conformal models has 
been given 0 § I 0 in terms of specific reductions of the two only existing two-dimensional single-boson 
integrable models, namely the sinh-Gordon and the Bullough-Dodd (BD) model (the Affine Toda Field The¬ 
ories and cSp) in their complex coupling constant versions which are also referred to as the sine-Gordon 
model and the Zhiber-Mihailov-Shabat (ZMS) model respectively SIMi- 

it is widely believed that the knowledge of the scattering data amounts to an exhaustive solution of a 
quantum held theory and that, in principle, one should be able to recover from them the operator content of 
the theory as well as to compute the correlation functions of local operators. In order to carry out this task, 
the so-called form factors approach has been developed and successfully employed in many important cases 
14 [l^, [Ig|, [IT], [is] Bf |2(], ||l], H|]. The strategy of this technique relies in the reconstruction of correlation 


functions by means of a spectral suni 1 ] 


1*1 (*) $ 2 ( 0 ) 10 >=£ 


dd i 
2tt 


M 

2tt 


n F*'(6u...,9 n ) [F**(6 1 ,...,6 n )]* e - m WZ* cosh<> ‘ 


( 1 . 1 ) 


on all the intermediate n-particle states of a scattering theory involving on-shell amplitudes of local operators 
(the so-called form factors) 

FfWi, ■ ■ ■ A) = (O|$(O)|A(0i)... A(0 n )). 


Form factors in turn can be exactly obtained in two dimensional bootstrap systems as solutions of a system of 
functional equations which entail the correct analiticity and monodromy properties dictated by the 5-matrix. 
The space of solutions of the system is supposed to give a faithful representation of the operatorial content 
of the theory j23| , but the correct identification of the form factors of a specific operator within this space of 
solutions is in general a nontrivial problem. A useful criterion for this identification was given in |is| ] where it 
was proved that for a scaling operator <F of scaling dimensions 2A$ the form factors divergence for large values 
of the rapidities is bounded by 

lim F*(0 1 ,...,0 n )<e A *l fl *l. 

\ 6 i |—>00 

■^For notational convenience we assume here the spectrum to consist of a single particle A and parameterize the momenta in 
terms of the rapidity variable 6. 
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More recently, it has been observed in |24|] that the form factors of relevant (A$ < 1 ) scaling operators satisfy 
a simple factorization property given by the so-called “cluster equations” 


A 


lim F®(6\ + A,..., 0 m + A, 0 m+ i,..., 9 n ) — F®(0i ,..., 0 m ) F®_ rn (9 m + 1 , • • ■, 0 n ), (1-2) 


F q 

1 n 


(Vm = 1 ,,n — 1), which hold unless some internal symmetry of the theory makes some of the form factors 
vanish. This property had already been noticed to be satisfied in the solutions of some specific models jj], 
[i7i[ p5| , pd[| and is believed to be a distinguishing property of exponential operators in Lagrangian theories. In 
particular, in ref. [ p5| a family of cluster solutions was found in the sinh-Gordon model and further identified 
with the form factors of exponential operators. These solutions were then used to compute the form factors of 
primary operators p7| | in a class of (f>i ^-deformed minimal models in which the boson of the original Lagrangian 
theory is still present after reduction. 

In the present paper we analyze the form factors of scalar operators in the Bullough- Dodd model and focus 


in particular on possible solutions of cluster equations (1.2). In view of the above discussion, these solutions 
become particularly interesting, not only as candidate solutions for the identification of exponential operators 
of the model, but also because in the complex coupling constant version (ZMS model) one should be able to 
identify among them the scaling relevant primary fields in the reductions which describe specific deformations 
of minimal models. 

The paper will be organized as follows. In Section 2 we review some general features of the BD model 
and its intepretation as a Complex Liouville Theory that allows us to map the exponential operators into the 
scaling primary fields of the <^ 1 , 2 , <t> 2,1 and (f>i^ deformed minimal models. In Section 3 we analyze the general 
solution of form factors equations for non-derivative scalar operators of the BD model, exhibiting the first 
multiparticle form factors and giving a full characterization of the dimensionality of the space of solutions. In 
Section 4 we study a one parameter family of cluster solutions which is shown to correspond to the form factors 


of exponential operators e kgv ^ and we determine the exact formula, eq. (4.10), which gives the dependence 
on k and g of these solutions. In Section 5 we check the validity of the correct interpretation of cluster solutions 
by comparing our results with all the known cases of form factors of scaling primary operators computed in 
(j> 1,2 and (j) 1,5 deformations of minimal models. In Section 6 we make use of the knowledge of the form factors of 
exponential operators to exactly compute the wave function renormalization constant of the BD model and the 
form factors of the fields <p(x) and :tp 2 (x): correctly normalized. Finally, we draw our conclusions in Section 7. 


2 The Bullough—Dodd model and its interpretation as a Complex 
Liouville Theory 

The so-called Bullough-Dodd (BD) model |jll| is a two-dimensional integrable Lagrangian QFT, namely the 
a^ 2 "* Affine Toda Field Theory Jl3| ], defined by the Lagrangian density 

£ = \ (« 2 - ^ (2 : e ^: + : e ~ 2 ™:) , ( 2 . 1 ) 
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where the exponentials are normal ordered. This model is the only 2D integrable QFT involving a single 
bosonic field which exhibits the p 3 property (i.e. the elementary particle appears as a bound state of itself). 
The spectrum of the theory consists of a single bosonic massive particle A of mass m and the exact R-matrix 
of the model factorises into two-particle amplitudes given by the following function of the relative rapidity 
variable 9 |29|| 


S(9)= / § (9)fs-,(9)f_B(d), 


where we have used the building block function 

_ tanh ^(9 + in x) 
x tanh i (9 — m x) 

The renormalized coupling constant B is given by 


B(g) 


g 2 / 2 7T 
1 + 3 2 /47t ’ 


and ranges from 0 to 2 for real values of g. The R-matrix exhibits a weak-strong coupling duality under 
the transformation g <—> 47 T /g or equivalently B <—> 2 — B. For later use we also define the following 
duality-invariant function of the coupling constant 


{B + 2) 7T 
c = cos--- 


( 2 . 2 ) 


The R-matrix has a simple pole at 9 = 2tti/3 corresponding to the bound state represented by the particle A 
itself. The on-shell three-point coupling constant is given by 

F 2 = -i lim (9-^)S(9) = 

0 -» 2 w /3 3 ' y (c-l)(l-2c)’ 

and vanishes both at the free field limiting values B = 0, 2 and at the self-dual point B = 1. 

For imaginary values of the coupling constant g (i.e. B < 0), the BD model — which is then referred to 
as the Zhiber Mihailov-Shabat (ZMS) model [|l2| permits the description of (f>ip and <p 2 ,i deformations of 
conformal minimal models. Indeed, starting from the observation that the ZMS has a non-unitary S'-matrix 
related to the Izergin-Korepin R-matrix, Smirnov exploited the quantum group SL{2) q invariance of the S'- 
matrix in order to recover unitarity in specific reductions of the model. The S-matrices of the above-mentioned 
deformed minimal models were in this way obtained from RSOS restrictions of the Izergin-Korepin R-matrix 
at specific values of the coupling constant at which q is a root of unity |lo|]. More recently the possibility of 
studying some relevant (j )\ } 5 deformations of specific non-unitary minimal models has been considered as well 
which relies again on quantum group reductions of the ZMS model [^0| . 

Some insight can be obtained if one considers the ZMS model by interpreting one of the exponential operators 


in the Lagrangian as a deformation of a Complex Liouville Theory (CLT) 131 . The other exponential then 
plays the role of a screening operator. If we require the CLT to describe the minimal model A4 r s with central 
charge 

6 (r - s) 2 

G = 1-. s > r relative primes , 

r s 
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and primary fields 0 TO)Jl of conformal dimensions 


A m,n — 


(nr — m s ) 2 — (r — s ) 2 
4r s 


m = l,...r— 1 ; n = l,...s— 1 , 


the above interpretation leads to a four-fold choice: in fact, after choosing one of the two exponentials in eq. 
dj as the screening operator, one can still choose two possible values of g as a function of r and s in order to 
correctly set its conformal dimensions to be A = 1. In the Complex Liouville Theory, the primary operators 
will be given by the following exponential operators'^ 


A — p k m, n g<P 


m= — 1; 77 = 1, ...s — 1; 


(2.3) 


with the identification 


u m,n — 'Vr—m,s—n ? 


(2.4) 


which entails the correct symmetry of the Kac Table of the model. In eq. (2.3), the dependence of fc m . n and 
B(g) from the integers r, s, to, n depends on the choice made for the screening operator, as summarized in 
Table 1 and the deforming exponential can be easily shown to correspond to one of the primary fields 01 , 2 , 02 ,l, 
or 05 .i. However, one can easily check that, while the primary field 0 1.2 is relevant in any minimal model, 
the field 0 5.1 is on the contrary always irrelevant and therefore does not yield renormalizable deformations. As 
for the fields <$> 2,1 and <^ 1 . 5 , they can be shown to be relevant only in disjoint sets of models: the field 02,1 is 
relevant for the class of minimal models M r .s with s < 2r which includes all the unitary cases A4 rir +i, while 
01.5 is relevant for the complementary class of non-unitary models s > 2r. Notice that in order to decide 
whether the deformation is relevant or not it is sufficient to require that the coupling constant g be imaginary, 
namely that B < 0 (see Table 1). 

The spectrum of the reduced ZMS model in general consists of kinks together with a cascade of their possible 
bound states [[To| . This spectrum does not always contain the original BD boson: while this particle is always 
present in the 01.2 deformations (where it appears as the lightest breather of two fundamental kinks), it is on 
the contrary never present in the spectrum of 02,1 deformations. In the relevant 0 i .5 deformations the presence 
of the BD boson depends on the specific model (see ]3C|). 

We will not enter in further details on the reductions of the ZMS model which can be found in the original 
literature ju], [io| . The information collected in this chapter is all we need for establishing the correct mapping 
between exponential operators of the BD model and primary fields of the reduced models. 


3 Form Factor Equations in the BD model 

We now turn to the problem of determining the on-shell matrix elements (form factors) of a local operator 
<h(:r) in the BD model. In the framework of two-dimensional integrable QFT, this problem is reduced to the 

2 In the following we will always consider normal ordered exponential operators omitting the notation : e a r ( x ) : . 
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F n (0 1 , • ■ ■ ,O i: 9 i+ 1, ...,9 n ) = S(6i - 9 i+ 1) F n (Oi,.. . ,0i+i,0i, ...,O n ), 
F n (d 1 + 2t ri, 0 2 • • •, 0„) = K(0 2 ■ • ■, 0 n , 0i), 


problem of studying a set of coupled functional equations 114, [15| namely the Watson monodromy equations 

(3.1) 

as well as the recursive residue equations on annihilation poles (kinematical residue equations) 

lim F n+2 (9' + iTT,9,d 1 ,...,0 n ) = i [l-]} 5(0-00) F n (0 1 ,...,0 n ) 1 (3.2) 


and bound state poles (dynamical residue equations) 

f 27 TZ \ 

Jirn^ [a-—) F n +2 (0 + a/2, 0- a/2, 0 ,,..., 0„) =iTF n+1 (Mi,---, 0„). (3.3) 

The most general solution to the monodromy equations (H|) can be written in the following form 0 

F*(6u. .., e n ) = R*{6 r,..., 0„) n F min (0* - 00 , 

i<j 

where R®(9 1, ..., 9 n ) is any symmetric 27ri-periodic function in the variables 9i and the “minimal” two-particle 
form factor F min {9) is given by the following function 

So (0) 03.(0) 


where g a (0) is defined by 


g a (9) = exp 


F mm (9)=Af(B)- 

g2-B_{9) gB{9) 

dt cosh ((a — 1/2)0 • 2/r 


(3.4) 


/o t cosh(t/2) sinht 


sin“((j7r — 9) t/2n) 


In eq. (3.4), A f{B) is the following normalization constant 

dt sinh(f/2) cosh(t/6) 


Af(B) = exp 


* 4 . 7 


(cosh(t/3) — cosh((i? — 1) t/ 3)) 


(3.5) 


sinh t 

chosen such that F mm (oo) = 1. For real values of the coupling constant, namely for B £ (0,2), F min (9 ) has 
neither poles nor zeros in the physical strip Im 9 £ (0,7r), since the same property is shared by g a (9) when 
a £ (0,1). The analitical continuation of F mm {9 ) for imaginary values of the coupling constant g (B < 0) 
developes poles in 9 which can be explicitly exhibited by using the following functional relations 


Si+a(0) = 9—a (0), 

g a {9)g- a {9)=V a {9) = 


cos 7ra — cosh 9 
2 cos 2 


satisfied by the functions g a {9). 

Notice that we have not mentioned yet the dependence of the form factors iqf on the operator $(x). 
Indeed, in the system of equations © and ( |3.3| ) this dependence is not explicit and further physical 

requirements are necessary to identify in the space of solutions the form factors of a specific operator. 
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3.1 General solutions for scalar non—derivative operators 

In this work we are mainly concerned with the analysis of form factors of scalar operators which are local 
nonderivative functions of the field tp(x). This infinite dimensional operatorial space can be spanned for 
instance by the basis of polynomials in (p(x) or by the basis of exponentials e av ^ x K A suitable parameterization 
of the form factors for this class of operators is the following 


rpmin f n . _ n \ 

F*(e u . ..,9 n )=H% Qt (Xl, ..., x n ) n ( Xi + x ) {x 2 + + x 2) ’ 


(3.6) 


i<j 


where Xi = e 6i . The pole structure expected to reflect the correct analiticity properties is explicitly shown in 
the denominator of ( |3.6| ), where annihilation and bound state simple poles are present at relative rapidities 
9ij = in and Oij = 2-7tz/ 3, respectively. is a homogeneous symmetrical polynomial in the variables Xi whose 
total degree is determined by Lorentz invariance to be d n = 3 n l" -1 ) . The constants H® are conveniently 
chosen to be 

H*=t^ n {B), (3.7) 

in order to obtain a simplified version of recursive equations on the polynomials Q In eq. (3d), t is a free 
parameter which will have an important role in the discussion of cluster solutions whereas 

V3T(B) 


KB) = 


J?min ^27rz ^ ' 


With the above choice of H ®, the dynamical recursive equations ( |3.3| ) read 

Qn(ux, u~ x x, Xi,..., x n - 2 ) = - x 3 D n _ 2 {x\xi,..., x n - 2 ) Q n -i(x, Xi,..., x n - 2 ), 
where ui = ! 3 and the polynomial D n is given by 


(3.8) 


D n (x\xi ,..., x n ) = ^ x 3n kl k2 k3 a ( ^ ] cos ((fc 2 - k 3 )(B + 2)tt/3) . (3.9) 

ki,k 2 ,k 3 =0 

The last expression is written in the usual basis of symmetrical polynomials cq n) which are defined by the 
generating function 

n n 

yi xn ~ k °i n)= n (*+ x *) ■ ( 3j °) 


fc=0 


i=1 


In expression (3.9) we can get rid of the trigonometrical dependence on the coupling constant B by exploiting 
the following recursive relation 

cos((n + l)a) = 2 cos(na) cos a — cos((n — l)a) , 

which allows us to express cosines of multiple angles as polynomials of cos a. In this way we can cast the 
dependence of eq. ([0]) on the coupling constant into a rational dependence on the variable c defined in eq. 

s 


The kinematical residue equations on annihilation poles (3.2) can be written as 

Qn( Xj Xj X \ , . . . , X n — 2 ) — ( ) K X U n — 2 {x | X\ , . . . , X n — 2 ) Qn— 2 (*£ 1 , - • - , X n — 2 ) , ( 3 . 11 ) 
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with 


U n (x\x u ...,x n ) = 2 ^2 (-) k2+k3+k *x 6n -( k '+'-- +k °K£ ) ---^ n 


k 6 

ki,...,ke=0 

• sin ((2 (k 2 + — ki — fc 3 ) + B (fc 3 + k e — fc 4 — k 5 )) tt/3) 


(3.12) 


and 


K = 


(2c-1) 


4 v / 3(l + c)(2c+l) ' 

Before solving the system of recursive equations, let us derive some important properties on the space of 


solutions from a direct analysis of the equations (3.8) and (3.11). 

A — It is easy to prove that in the space of symmetrical polynomials of degree d n = 3 n (” . ' ] ) ) the only 
polynomials which have zeros both at Xi/xj = e 2?ri / 3 and at Xi/xj = —1 are given by 

K,^ n \{xi}) = (xi +xj)(xj +XiXj +Xj) 

l<i<j<n 


= det 


T («) 
7 2 j-i 


det 


.(«) 

z[j/2\-i+i+{-Y 


l<i,j<2n-2 


up to a multiplicative constant. This is therefore the only possible kernel for the whole system of recursive 
equations. Hence, after fixing all the polynomials Qi for i = 1,... n — 1, the most general solution Q n of 
the system of equations (3.8) and ( |3.11 ) will be then given by 


Qn = Q* n +XnlC in \{Xi}) 


(3.13) 


where Q* n is a specific solution and X n is a free parameter. The space of solutions will be organized 
correspondingly, namely every operator will be identified by a succession of parameters \, i = 1,.. .oo 
and the general solution for a n-particle form factor will be described by an n-dimensional vector space 
of solutions Q n spanned by the parameters Ai,..., A„. 

(i) 

B — The partial degree of the general polynomial Q n with respect to any of the variables ais exactly d n = 
3(n — 1). This can be easily shown by induction observing that Q\ must be a constant for Lorentz 
invariance and making use of equations (3.8), (3.11) and ( 3.13| ). This implies in particular that the form 
factors of this class of scalar operators of the theory have bounded asymptotic behavior for large values 
of the rapidities, 


lim F*( 6 1 + A,..., 6 k + A, 6 k + 1 ,..., 0„) < oo 

A—K50 


Vfc = 1,... ,n — 1. 


This observation enables us to look for cluster solutions of form factors equations within this general class 
of solutions (see eq. dU))- 

We now turn to the actual computation of the first multiparticle general solutions to the system of recursive 


equations ([hSl) and ( 3.11 ). The most direct way of computing these solutions consists in parameterizing any 
polynomial Q n as the most general polynomial of degree d n = 3 ” l” -1 ) i n the basis of symmetrical polynomials 
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aj' n> and to impose on the coefficients of the expansion the costraints coming from the recursive equations. 
We report here the result of the first general multiparticle form factors in the space of scalar non-derivative 
operators. Lorentz invariance requires Q± to be a constant 

Qi = Ai , 


hence in order not to have two different overall normalization constants we can set for the time being t = 1 in 


eq. (3/7). The next most general solutions are given by 

<92(2:1, *2) = -Aicri 3 - A 2 /C (2) , 


<93(^1, £2,3:3) = 


. / 4 4 ( 4 C 2 - l) 2 2 

Ai CTi cr 2 + C1 cr-2 tr 3 + ——-- ai 02 CT3 — 


( cr 2 3 tr 3 + 01 3 cr 3 2 ) 


2 (1 + c) * " 0 2 (1 + c) 

+A 2 (<t 1 cr 2 4 + CTi 4 ct 2 ct 3 — 2(1 — c) ( 01 2 ct 2 2 ct 3 — 01 cr 2 <r 3 2 ) - (J 2 3 cr 3 - di 3 cr 3 2 ) 

+a 3 /c( 3 >, 


where the residual kernel freedom of each solution has been explicitedQ. Notice that in the above solutions the 
trigonometrical dependence on the coupling constant has been hidden in a simple rational dependence on the 
self-dual variable c defined in eq. (|2.2|). This major simplification has been made possible by noticing that the 


systematic solution of the dynamical recursive equations alone (3.8) yields polynomials Q n which already have 
the correct single-parameter kernel ambiguity ( 3.13| ) expected for the whole system. It therefore means that, 
actually the dynamical recursive equations (3.8) are equivalent to the system of the two coupled equations 


and (3.11). 

The general solutions that we have found must include in particular the form factors of the elementary field 
<f(x) which were first studied in p8[| . One can prove that they can in fact be selected by imposing either the 
asymptotic vanishing of the form factors for large values of the rapidities (i.e. imposing the cancellation of the 
highest partial degree terms in Q n ) or the proportionality Q n ~ a n j32|. The A i are then determined to be in 
this case 

\V _ _ \V 

^2 — , 


Af = 0 V*> 2. 

Finally the overall normalization is fixed by|j] 

Z 1 ! 2 

< 0 | v (oM = -7p 

which sets A) 3 = /r _1 Z 1 / 2 /v / 2. In the above expression Z is the wave function renormalization constant of the 
theory which will be exactly computed in Section 6. 

3 We do not report here the general solution of Q 4 which already contains an extremely large number of terms and is not 
particularly useful for the purposes of this work. 

4 Our convention on the normalization of states is (A{0i)\A(02)) = 2ttS(0i — 62) = 2ttE\ 5(pi — P 2 ) • 
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By using the above general solutions we can also identify the 1-parameter family of the trace 0(:r) of the 
Stress-Energy tensor for different values of the background charge. This family of operators was studied in ref. 
pfil where the authors showed that different choices of Q(x) select different possible ultraviolet limits of the 
theory. In order to identify these form factors it is sufficient to impose the proportionality Q n ~ <j\ a n -\ for 
n > 3, as it can be shown from the conservation of the Stress-Energy-Tensor. In this way one determines all 
the free kernel parameters A j but the first two. The parameter Af is found to be for example 


A 


e 

3 



3 Af 
2 c ~h 2 


Finally, imposing the overall normalization 


F ’® (in) = 2 7r to 2 , 


one determines 

x©_ Km 2 

2 (c -1) r 2 ’ 

and obtains a one-parameter family of independent operators for arbitrary A® which coincides with the one 
analyzed in ref. p(i|| . 

In order to identify different operators in this general space of solutions one must resort to more powerful 
techniques. We will see in the following section how the imposition of the cluster equations (jb^) enables us to 
extract the form factors of a whole basis in the space of non-derivative scalar operators. 


4 Form Factors of Exponential Operators 


In this chapter we study the existence of solutions of the form factor equations which also satisfy the further 
requirement given by the so-called cluster equations fll.2| ) imposed on a multiparticle form factor F n . This 
restrictive set of non-linear equations is believed to select out the exponential operators in a Lagrangian theory 


15, ^5|. More recently it has been shown in ref. [g4| that these equations are the distinguishing property of 
scaling operators in the conformal limit of a two-dimensional field theory at least in the cases where there is 
no symmetry preventing the form factors from being non-vanishing. This observation has been confirmed and 
successfully employed in ref. ’s §1 for identifying the complete set of scaling primary fields in some massive 
deformations of minimal models. Cluster solutions become therefore objects of utmost interest in the BD 
model because the two ways of looking at them either as exponential operators or as scaling fields, converge in 
this theory where specific exponentials are identified with primary operators in the reduced models describing 
deformations of conformal field theories. 

In order to impose the cluster equations (0 we fix the overall normalization of the form factors by adopting 
the convenient choice F 0 = 1 and choose 

Qi = 1. 
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Equations (|l^) then amount to requiring the following property on the polynomial Q n 


lim 


Qn 


Qn 


Qn—r. 


A—too K.( n ) ^ 2 ' 1 ’ ' ' " ’ ^ Xm ' 1 ■ ■ ■ > X n) — ^ . ,,ra) 


(4.1) 


where t - the variable introduced in eq. & - is now switched on and treated as a free parameter. These 
further restrictions imposed on the general solutions of residue equations determine level by level all the X n 
parameters as functions of t. At any given level n , the number of equations which determine the only free 
parameter left X n , grows rapidly with n, therefore the very existence of a cluster solution is not at all obvious. 
For the first computed solutions however, all the equations on a given X n turn out to be identical and we 
believe that this should be the case at any level. In this way we obtain a one-parameter family of solutions for 
t arbitrary, of which we report the first multiparticle representatives]^ in Appendix A. Notice that t is not an 
overall normalization factor since the normalization of the form factors has been fixed by Fq = 1 and indeed, 
due to the nonlinearity of (0, the solutions Q n (t) turn out to be polynomials in t of degree n — 1. This 
means that t defines through the polynomials of Appendix A and eqs. 


and (3.7) a one-parameter family 
of solutions Fn ^ corresponding to independent operators. If we make then the hypothesis that these solutions 
actually correspond to the form factors of the exponential operators e kav ^ x \ 


p{ t} m a , _ (0\e ka ^\A(6 1 )---A(6 n )) 
n (0|e fc ®^(°)|0) 


(4.2) 


we are forced to consider t as a well-defined function t(k,B) of k and B rather than a free parameter. In 
particular, in order to establish the one-to-one correspondence between cluster solutions and exponential 
operators it is of particular interest to compute the normalization-invariant quantity 

'"\ e kg<p(0)\ A } 


p 1*1 = 

1 (0|e fe ^(°)|0) 


= v{B)t{k,B) . 


(4.3) 


In the following we consider in detail some conditions that we can impose on the function t(k,B) in order to 
find its exact form. 


4.1 The Function t(k, B) 

The first information on t(fc, B) can be obtained from the computation of the conformal dimensions 
A = — g 2 k 2 /Sty of the operators e kav ^ x ' > in the free-boson ultraviolet limit at lowest order in g 2 . These 
can be easily obtained from the analysis of the short distance behavior of the correlator (0|e fc9¥ ’^' ) e fc9V ^ O ' ) |0) by 
means of eq. 0 ) and the cluster solutions F^ ' 1 . We obtain 

A = - 2 iim = 9 2 t(k, 0) 2 

^ g—>o 4 7 t g 2 8 ir 

from which one obtains the important relation 


lim t(k. B) = k . 
B —>0 


(4.4) 


5 Notice that from a computational point of view there is no difficulty in obtaining the next multiparticle solutions since the 
dynamical recursive equations (3.8) are linear equations in the unknown coefficients of independent monomials in the cr’s and the 


dependence on the coupling constant is simply a rational dependence on c. 
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Furthermore, from the expressions (A.S) and (A.4), by imposing the proportionality Q n ~ <Ji<r n -i, one can 
easily verify that the only cluster solutions which also belong to the class of possible traces of the stress-energy 
tensor are defined by the solutions of 


-l + 2c + 2t + 2ct + 2t 2 + 2ct 2 = 0, 


namely 


= 


sin((i? + 1)7 t/ 6) 
cos ((B + 2)7t/6) 

sin ((B — 3)7t/6) 


(4.5) 


cos {{B + 2)7t/6) 

These two solutions correspond to the ones found in ref. |2tJ and identified with the form factors of the 
fundamental vertex operators e 9V and e~ 2gv which appear in the Lagrangian density. This can also be obtained 
immediately by taking the limit B —> 0 in eq. (4.5) which gives respectively k = 1,-2 in virtue of (4.4). 
Therefore we have also the two following important requirements on t(k, B): 

sin((.B + 1)7t/6) 


t(-2,B) = 


cos ((B + 2)7t/ 6) ’ 
sin ((£? — 3)7r/6) 


(4.6) 


(4.7) 


cos ((B + 2)7t/6) 

As a limiting case of the cluster solutions we can also recover the form factors of the fundamental field <p(x) 
which is naturally obtained from the vertex operators in the limit k — > 0. These form factors of course satisfy 
a trivial cluster property because they vanish for large rapidities and therefore satisfy eq. Q with t = 0. 
Hence we get one more information 

lim t(k, B) = 0 . (4.8) 

k —>0 

Indeed one can easily check that the form factors we had obtained in the previous section for the held ip{x) 
from the most general solutions of residue equations satisfy, 


F% = Af lim 

,L 1 4- , n 


F, 


{*} 


i—*0 t 

A remarkable check on the correct identification of these operators is obtained studying the quantum equations 
of motion of the model 


m n 


Utp + ^ ( e 9V - e~ 2av ) = 0. 


3 9 


If our identihcation is correct we should hnd[] 


m 


T («) 
T n-1 


(n 
CTn 


FX + t = 0 , 


with some constant r, or equivalently 

(n) (n) 

Afm 2 1 




' Qn(0) + T (t+ Q n {t + ) - t~ Qn{t~)) = 0 . 


6 In general J 7 ')) 14 ’ = —m 2 <71<T " 1 F* for any field <£(2:). 
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Indeed this last equation can be verified to hold on the solutions given in Appendix A with 


\ v 2 

t =-^=- tan((B + 2)tt/6). 

The non-perturbative nature of this last check shows that the identification of cluster solutions as vertex 
operators is far beyond a semiclassical one for small coupling constant. 

The constraints obtained for the function t(k,B), eqs. (4.4), ©, (O) and ( [4.8[ ) are not sufficient to 
determine its form and, in particular, little information is given on the dependence on k. We will see however 
in the next section that some additional requirements coming from the reductions of the ZMS model impose a 
periodicity condition in k for the function t(k,B) 


t{k , B) = t{k + 6/13, B ), 


(4.9) 


which suggests the following conjecture: 


t(k, B ) 


sin(fcI3 7r/6) sin((k B + B + 2)n/6) 

2 sin (B 7t/ 6) sin((2 — B)tt/6) cos((i3 + 2)7 t/6) ’ 


(4.10) 


This function satisfy all the aforementioned requirements. A decisive check of the validity of this expression will 
be obtained in the following chapter by the comparison with explicit computations of form factors of primary 
operators in specific reductions of the ZMS model. This formula may be regarded as one of the main results 
of the paper. In fact it allows us to explicitly assign to every vertex operator e kav in the BD model its form 
factors Fn k ^ which are obtained from the cluster solutions Q n (t) of Appendix A through the parameterization 
( pd^ ) and eq. (0) by replacing t = t(k, B ). 


5 Form Factors in the reductions of the ZMS model 


We now turn our attention to the analitical continuation of the model to imaginary values of the coupling 
constant g, namely to possible reductions of the ZMS model. In these models the spectrum is no more a single¬ 
particle one as in the real coupling BD model, but it has a richer structure that depends on the model analyzed. 
We consider here only those restrictions whose spectrum still contains the elementary boson excitation of the BD 
model, namely 0i,2 and some cases of </p 5 deformation^. If we assume that the identification obtained between 
cluster solutions and vertex operators of the model is exact, we are then led to establish a correspondence 
between the form factors of exponential operators e k9V ^ in the BD model and the form factors^ of scaling 
primary operators in the deformations according to the correspondence given by eq. (|2.3| ) and Table 1. An 
immediate consistency requirement for this procedure is obtained by imposing that the form factors respect 
the symmetry (2.4) of the Kac table of minimal models. For example, the quantity of eq. (4.3) should 


have the same value if evaluated at k = fc m ,n and k = fc r _ mjS _ n . Imposing this condition both in the <f>i^ 

7 To avoid confusion we stress that the elementary BD scalar boson, which is created from the field ip(x), is not the fundamental 
particle in the bootstrap of the reductions of the ZMS model which is instead a three-component kink. 

8 The form factors in the reduced model must be intended as the matrix elements on the BD breather sector. 
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deformations and in the (p\$ relevant ones we obtain respectively that the following two symmetries of the 
function t(k,B) must hold 


t(k, B) = t{-k-l-2/B,B), 
t(k, B) = t(—k — 1 + 4/13, B ), 


(5.1) 


which in particular entail the above mentioned periodicity in k, equation (4.9). Both these symmetries are 


indeed separately satisfied by the function ( [4.101) . 

A precise check on the validity of equation ( 1.10| ) is provided by comparing its predictions with the form 
factors of scaling primary operators in (pi ^ and deformations which can be found in literature. We have 


indeed computed the normalization invariant ratio F\/F$ using eq. (4.3) and the assignments of Table 1, for 
all the known cases of primary form factors which have been analyzed in literature |LG|, |l9[ ^l| [20| Q (see 
Table 2) and a perfect agreement has been found with all the values reported in the references. We stress 
here the fact that in the references considered, the form factors of primary operators have been identified by 
different techniques: in ref.’s [|I(| |F], the identification has been obtained by using the correspondence 
between the deforming field and the trace of the stress-energy tensor whereas in ref. ’S01 the form factors 
of the primary fields have been identified with the finite number of solutions of a non-linear system of cluster 
equations involving the form factors relative to the whole particle spectrum of the reduced models. 


6 The Wave Function Renormalization Constant of the BD Model 
and the Form Factors of cp(x ) and : ( p 2 {x): 

The form factors Fr) k ' > {6\ ,..., 9 n ) that we have computed have been so far conveniently normalized putting 
1*0 = 1. From equation ([b^) one immediately observes that these form factors are invariant under an additive 
redefinition of the field <p(x) — > ip(x ) + const. We remove this ambiguity on the definition of the field ip(x ) by 
imposing that its vacuum expectation value (0|<p(:r)|0) be zero, namely subtracting from the original Lagrangian 
field the value of the one point tadpole function. Consider now the following expansion of the form factors of 
exponential operators: 


^ 7,jf j 

<0|e fc ^(°)| A{6 1 ) ■ • • A(B n )) = ]T —^<0| :^(0): | A(0 1 ) ■ ■ ■ A(0 n )), 

j =1 

and of the vacuum expectation value 

' h,j j 

<0|e fc ^(°)|0> =£-r<°l w'(0): |0> = 1 + o(k 2 ). 
j =0 J - 
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If we now expand the form factors that we have obtainedjj in series of k we can identify the form factors 
of the fields p(x) and '-F 2 ( x ) : as the coefficients of order k and k 2 respectively. 

M 0 )\A{e 1 ) ■ ■ ■ A{e n )) 

( 6 . 1 ) 




, #n) = 


(0| e fc 9¥ .(°)|0) 


2 „2 


= kg (O|y)(a;)|A(0i) • • • A{9 n )) + 


k<J 


:tp 2 (x): \ A(6i ) • • • A(0 n )) + o(k 3 ). 


This procedure gives the form factors of <p(x) and :<p 2 (x ): with the correct overall normalization of the fields. 
This observation in particular allows the exact determination of the wave function renormalization constant 
Z(B) of the BD model. In fact, considering the first order expansion in k of 


F, 


f> = MB) i(k, B ) = „(fi) , „ B ,1 T' ((g APS: 


12 sin(i?7r/6) sin((2 — B)ir/6) 


,{k 2 ) 


= ^ ^ (0|^(0)|m4) +o(k 2 ) 


kgZ 1 / 2 


+ o(k 2 ). 


one easily obtains the following expression for Z(B) 
Z{B) = ^B) 2 B (2 - B) 


tan((i? + 2)7 t/6) 


288 \sin(i?7r/6) sin((2 — B)n/6) 
2 tt B(2-B) (c — 1 ) 


( 6 . 2 ) 


3 y/3 A f(B) (1 + 2 c) (1 — 2 c) 

where Af(B) is defined in eq. (|3.5| ). The function Z(B) is manifestly dual with respect to the weak-strong 
coupling transformation B <—> 2 — B and can be easily shown to coincide at lowest order in g 2 with the correct 
perturbative result coming from the one-loop self energy diagram 


Z = l- 


9 


1 


1 


°(g 4 


12 3V3y 

A plot of the function Z(B) is given in Figure 2. Notice the tiny deviation of the constant from the free field 
value Z = 1 on the entire range of the coupling constant B £ [0, 2]. 

The correctly normalized form factors of the field p(x) are given by 


F V - cr 1 — F W 
" “ 9 dk n 

Z x ' 2 F^ t} 


(6.3) 


k—0 


t 


t -o 


while the exact form factors of the field : ip 2 (x ): are simply obtained by 


For example we can compute 

f y = 


F A 2 = a -2 F {k} 
9 dk2 *n 


:p 2 (0):\A) 


(6.4) 


k =0 


9(B) 9 2 


fe=o 


= f{B)B(2-B) 


1 


144 sin(B7r/6) sin((2 — B)tt/6) 


9 Here and in the following we will adopt the notation f}/'' instead of F ,| t ^ to stress the dependence on k. The relation between 
the two expressions is obviously given by t = t(k, B) eq. (ll. 1 Of). 
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which exactly matches at lowest order in g with the one loop calculation 


In a similar way we get 


(0|:^ 2 (0):| A)=-^= + o(g 3 ). 


F^\d 1 -9 2 ) = (0| :^ 2 ( 0 ): \A{0{) A(0 2 )) 

2 7T 1 

^ ^ ^ ^ ^ 288 (sin(f? 7 r/ 6 ) sin((2 — f?)- 7 r/ 6)) 2 

• tan 2 {{B + 2)7 t/6 ) — a\a 2 {2 sin(I? 7 r/ 6 ) sin((2 — B)n/6) + tan 2 {{B + 2)7 t/6 ))^ 

F min ( 6 1 - 0 2 ) 

{x 2 + X 2 )(xj + XiX 2 + x%) ' 


Notice that in order to obtain the form factors of arbitrary operators :tp n (x ): one should exactly compute the 
vacuum expectation value {0\e kgv> ^ | 0 ) of the exponential operators and make use of expansion 0 ) (for the 
sine—Gordon model the vacuum expectation value of the exponential operators has been recently obtained in 

ref. JH ). 


7 Conclusions 


In this paper we have computed, in the framework of the bootstrap approach to integrable models, the first 
multiparticle solutions of form factor equations for general non-derivative scalar operators in the BD model. 
Among these solutions we have selected a one-parameter family of cluster solutions which have been identified 


by means of the central result eq. (4.1C) with the form factors of exponential operators e kav> . In the complex 
coupling constant version of the model, the form factors of exponential operators allow to identify the form 
factors of relevant primary operators in the sector of the lightest breather of <t>\_ 2 and cf> 15 deformations of 
minimal models and perfect agreement has been found with all the examples that we have found in literature. 
Finally, by using the cluster solutions, we have computed the form factors of the fields tp(x) and : <f 2 (x) : 
with the correct overall normalization and determined in this way the non-perturbative exact wave function 
renormalization constant of the model. 

We have therefore obtained the characterization of form factors for a whole basis in the space of scalar 
non-derivative operators and we have found complete consistency, in a non-perturbative setting, between the 
axiomatic ^-matrix approach to bootstrap systems and the Lagrangian approach to quantum field theories. 

This work also yields an efficient tool for the identification of relevant primary fields among the cluster 
solutions of massive <j>\ t2 and <j> 1,5 deformations of minimal models. 


15 






Acknowledgments: I wish to express my grateful thanks to G. Mussardo for his continuous and warm 
support in the preparation of this paper. I want to thank also G. Delfino, A. De Martino and F. Smirnov 
for precious discussions. I am grateful to the INFN for supporting a monthly stay in Paris at Institute Henry 
Poincare where the work has been completed and to prof. O. Babelon for his ospitality. 


16 



Appendix A 


In this Appendix we list the first solutions of the one-parameter family of Q n polynomials of cluster solutions 
in the BD model. In the following expressions, the variable c is the dual-invariant function of the coupling 


constant defined in eq. (2.2) and t is a free parameter. The solutions are identified with those of the basis of 
operators e k9V by means of eq. (4. 10|) which determines t as a function of k and g. 




(A.l) 


Q 2 (t) = ta\ (A.2) 

- (1 + t) 01 02 , 

Qz{t) 2 (1 + c) = 2 (1 + c) t 2 0 i 3 02 3 (A.3) 

— 2 (1 + c) t (1 + t) <7\ 02 ^ 

-2 (1 + c) t (1 + t) 01 4 02 03 
+ (3 + At — 4c 2 t — 2t 2 — 2cf 2 ) 0 i 2 02 2 03 
+ (— 1 + 2c + 2 1 + 2 ct + 2 1 2 + 2 ct 2 ) 02 3 03 
+ (-l + 2c+2f+ 2cf + 2f 2 + 2ct 2 ) 0i 3 0 3 2 
+4 (-1 + c) (1 + c) (1 + t) 01 02 03 2 , 

Qi(t) 2 (1 + c) = (A.4) 

= 2 (1 + c) t 3 01 3 02 3 03 3 

— 2 (1 + c) f 2 (1 T t) 0102 ^ 03 ^ 

— 2 (1 + c) t 2 (1 + f) 01 ^ 0203 ^ 

+t (3 + At — 4c 2 t — 2t 2 — 2 ct 2 ) 0 i 2 02 2 03 4 

+f (—1 + 2 c + 2 f + 2 ct + 2 t“ + 2ct 2 ) 02 3 03 4 

+f (—1 + 2 c + 2 t + 2 ct + 2 t“ + 2 ct ) 01 03 3 

+4 (-1 + c) (1 + c) t (1 + t) 0 i 0 2 03 5 

-2 (1 + c) f 2 (1 + t) 01 3 02 4 03 04 

+ 2 (1 + c) t (1 + t) 2 01 0 2 5 03 04 

+t (3 + At — 4c 2 t — 2t 2 — 2 ct 2 ) 0 i 4 02 2 03 2 04 

+2 (1 + t) (—2 + c — 2t + 2ct + Ac 2 t + 3t 2 + 3 c t 2 ) 0 i~“ 02 3 03 2 04 

+ (l + f) (\~2c~2t~2ct-2t 2 -2ct 2 ) 0 2 4 0 3 2 04 

+f (—1 + 2c + 2f + 2cf + 2t” + 2ct ) 0 i 03 3 04 


17 






+ (l — 2 c — 4i + 4 c 2 1 — 2t 2 + 14:ct 2 + 8 c 2 t 2 — 8 c 3 i 2 + 6 i 3 + 6 ci 3 ) c 4 3 c 2 c 3 3 c 4 

-h (7 — 8c -h 9 i — 14 ci — 12 c 2 i -h 8c 3 i — 2 i 2 — 6 ci 2 — 4 c 2 i 2 — 2 i 3 — 2 ci 3 ) c 4 C2 2 C3 3 <74 

-hi (3 — 14 c + 8 c 3 — 6 i — 14 c i -h 8 c 3 i — 6 i 2 — 6 c i 2 ) ci 2 03 4 c 4 

+2 (— 1 + c) (l — 2c — 2t — 2ct — 2t 2 — 2ct 2 ) 02 CT 3 4 c 4 

"hi (—l-h2c+2i + 2ci + 2i" + 2ci“) c 4 4 C2 3 c 4 ~ 

-h(l-ht) (l — 2c — 2t — 2ct — 2t 2 — 2ct 2 ) Ci 2 C2 4 a 4 2 
+4 (-1 + c) (1 + c) t (1 + t ) ci 5 c 2 c 3 c 4 2 

+ (7 — 8c + 9i — 14 ci — 12c 2 1 + 8c 3 t, — 2t 2 — 6 ct 2 — 4c 2 i 2 — 2i 3 — 2ci 3 ) c 4 3 c 2 2 c 3 a 2 
+ (—5 + 14 c — 8 c 2 — 6i — 2ci + 4 c 2 t — 6i 2 — 14ct 2 + 8c 3 1 2 — 4i 3 — 4ci 3 ) c 4 c 2 3 c 3 a 2 
-hi (3 — 14c + 8c 3 — 6i — 14ci -h 8c 3 i — 6i 2 — 6 ci 2 ) c 4 4 C3 2 a 2 
+2 (—5 + 14 c — 8 c 2 — 2t + 8 ct — 6 c 2 1 — 8 c 3 1 + 8 c 4 1 + t 2 — 11 ct 2 — 4 c 2 1 2 
+8 c 3 1 2 — 3 1 3 — 3 ct 3 ) c 4 2 c 2 C3 2 a 2 
+ (l + 4 c — 4 c“ + 2 i) (—l + 2c-h2i-|-2ci-h2i -{- 2 c i ) C2~ C3 c 4 

+ (—4 c -h 12 c 2 — 8 c 3 — 3 i -h 22 c i — 16 c 3 i + 6 i 2 + 22 c i 2 — 16 c 3 i 2 + 6 i 3 + 6 c i 3 ) c 4 C 3 3 c 4 2 
+2 (—1 + c) (l — 2c — 2t — 2ct — 2t 2 — 2 ct 2 ) ci 4 c 2 c 4 3 
+ (l + 4c — 4 c 2 + 2 i) (—l + 2c-h2i-h2ci-h2i 2 -h2ci 2 ) ci 2 c 2 2 c 4 3 

+ (—4c -h 12 c 2 — 8 c 3 — 3i -h 22 ct — 16 c 3 1 + 6 i 2 + 22ct 2 — 16c 3 i 2 + 6 i 3 + 6 ci 3 ) c 4 3 c 3 a 3 
+ (9 — 30c-h 20c 2 + 16c 3 — 16c 4 -h 8 i — 16ci+ 8 c 2 i + 16c 3 i — 16c 4 i 
+2 i 2 + 10 c i 2 — 8 c 3 i 2 + 2 i 3 + 2 c i 3 ) c 4 c 2 C 3 c 4 3 
+ (4c —4c 2 + i) (l — 2 c— 2 i — 2 ci— 2 i 2 — 2 ci 2 ) c 3 2 c 4 3 
-h(4c —4c 2 + i) (l — 2c — 2i — 2ci — 2i 2 — 2ci 2 ) c 4 2 c 4 4 . 
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Table Captions 


Table 1 Complex Liouville Theory assignments between exponential operators and primary fields for different 
choices of the screening operator. 

Table 2 Primary operators in ZMS reduced models for which the form factors have been computed in litera¬ 
ture. 


Screening operator 

Deformation 

B 

fcm,n 

e~ 2gv> 

e" 2 gv> 

e gip 

e 9 <P 

e 9<p = 0i i2 

e 9V = 02,1 

e~ 2g<p = 0 lj5 

e~ 2g(p = 05,i 

2 r 

r— 2 s 

2 s 

s— 2 r 

4 r 

2 r—s 

4 s 

2 s—r* 

(n- 1) — (to— l)f 

(to- 1) - (n- 1)| 
i((i-n)-( 1 -TO)^) 

3 ((! - m ) - l 1 ~ n )l) 


Table 1 


Model 

Deformation 

Primaries analyzed 

*1/^0 

Reference 

■M2,5 

01,2 

01,2 

0.8372182 z 

@ 

■M 2 ,7 

01,2 

01,2 

0.8129447* 

@ 



01,3 

1.245504* 

H 

M2 ,9 

01,2 

01,2 

0.7548302 * 

n 



01,3 

1.288576* 

n 



01,4 

1.564863* 

Hi 

M3,4 

01,2 

01,2 

-0.6409021 

©0 



02,1 

-3.706584 


•Ad 4,5 

01,2 

01,2 

-0.8113145 


Me,7 

01,2 

01,2 

-0.9499626 

@ 

M2,9 

01,4 = 01,5 

01,2 

-0.5483649 

n 



01,3 

-1.476188 

Hi 



01,4 

-2.169493 

H 


Table 2 
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Figure 1: Plot of the wave function renormalization constant Z(B) of the Bullough-Dodd model. 
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